Ab initio vibrations in nonequilibrium nanowires 
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Abstract. We review recent results on electronic and thermal transport in two different quasi 
one-dimensional systems: Silicon nanowires (SiNW) and atomic gold chains. For SiNW's we 
compute the ballistic electronic and thermal transport properties on equal footing, allowing us 
to make quantitative predictions for the thermoelectric properties, while for the atomic gold 
chains we evaluate microscopically the damping of the vibrations, due to the coupling of the 
chain atoms to the modes in the bulk contacts. Both approaches are based on the combination 
of density-functional theory, and nonequilibrium Green's functions. 



1. Introduction 

Modeling charge transport in nanostructures under nonequilibrium conditions via ab initio 
methods is nowadays a well-established computational approach [T]. Electronic transport, 
however is intimately coupled to the motion of the ions comprising the nanostructure: the 
charges may excite ionic motion which leads to local heating and energy relaxation. On the 
other hand, the ionic motion itself is an important channel for heat transport between the 
electrodes. The vibrations within the nanostructure couple to the motion of the atoms in the bulk 
contacts, thus becoming damped or, equivalently, they obtain a finite life-time. The inelastic 
transport properties of nanostructures is a very active research field, where progress occurs 
on many fronts. In this paper we review two particular subfields in this general area, (i) We 
compute the ballistic electronic and thermal transport properties of Si-nanowires, paying specific 
attention to intentional nanostructuring. The main goal is to examine whether the thermoelectric 
properties can be enhanced via nanoengineering. Treating electronic and vibrational properties 
on equal footing (based on nonequilibrium Green's functions) is a key element in our work. 
Presently it is not computationally feasible to consider inelastic scattering in these systems 
which typically consist of thousands of atoms. The important issue is thus to find a good 
compromise between microscopic accuracy, and the macroscopic properties of real, nanowire 
based devices, (ii) In smaller systems, such as the experimentally studied atomic gold wires, 
a fully microscopic approach has become tractable. The vibrational motion excited in the 
atomic chains due to electron transport is damped because the localized vibrations couple to 
the vibrations in the electrodes surrounding the gold chain. It is commonplace to describe this 



damping by phenomenological life-times; here, however, we describe a recently developed ab 
initio method to actually compute these life-times. Experimentally, a very rich behavior as 
a function of the elongation of the wire is observed, and the extreme sensitivity to stress, as 
revealed by our calculations, may form an important ingredient in explaining this behavior. 

2. Ab initio modeling of electronic and thermal transport in Silicon nanowires 
(SiNW) 

2.1. Background remarks 

Recent ground-breaking experiments indicate that rough silicon nanowires (SiNWs) can be 
efficient thermoelectric materials although bulk silicon is not [21 [3]: they conduct charge well 
but have a low heat conductivity. A measure for the performance is given by the figure of merit 
ZT = G e S 2 T/K, where G e , S, and T, are the electrical conductance, Seebeck coefficient, and 
(average) temperature, respectively. The heat conductivity has both electronic and phononic 
contributions, k = n e + K p h. Materials with ZT ~ 1 are regarded as good thermoelectrics, 
but ZT > 3 is required to compete with conventional refrigerators or generators [JJ. Recent 
theoretical works predict ZT > 3 in ultra-thin SiNWs [SEltll- The high performance SiNWs 
in Ref. [2] were deliberately produced with a very rough surface, and the high ZT is attributed 
to increased phonon-surface scattering which decreases the phonon heat conductivity, while the 
electrons are less affected by the surface roughness. The extraordinary low thermal conductivity 
measured in rough SiNWs is supported by recent calculations [HI [9]. Surface disorder will, 
however, begin to affect the electronic conductance significantly in very thin wires |10j and 
thereby reduce ZT [6j [7j . 

For defect-free wires, the room temperature phononic conductance scales with the cross- 
sectional area [11]. Contrary, for ultra-thin wires with diameters D < 5nm, the electronic 
conductance is not proportional to the area: it is quantized and given by the number of states 
at the band edges. Decreasing the diameter to this range thus decreases K p h while keeping 
G e almost constant. An optimally designed thermoelectric material would scatter phonons but 
leave the electronic conductance unaffected, even for the smallest wires. 

2.2. Optimizing the nanowire design 

Instead of introducing surface disorder, which affects the electronic conductance in the thin 
wires, we have recently addressed the question whether other surface designs could lead to 
improved thermoelectric performance [12] . In a similar spirit, Lee et al. [13J proposed already 
earlier nanoporous Si as an efficient thermoelectric matrial. Blase and Fernandez-Serra [14J have 
shown that covalent functionalization of SiNW surfaces with alkyl molecules leaves the electronic 
conductance unchanged, because the molecular states are well separated in energy from the 
nanowire bandedges, see Fig. [TJ If the alkyl molecules scatter the phonons, such functionalized 
SiNWs would be promising candidates for thermoelectric applications. Experimental alkyl 
functionalization of SiNWs was reported in Ref. [15j . 

Another possible surface design are branched SiNWs, so-called nanotrees. Nanotrees have 
been synthesized in III- VI semiconductors |16t [T7] and in silicon |18|. 119] . The stability and 
electronic structure of silicon nanotrees have recently been addressed theoretically [20[ [2TJ . The 
thinner branches will have a larger band gap than the main wire (" trunk" ) , see Fig. [TJ Close to 
the band edges, the electronic scattering is therefore weak. 

The presence of an alkyl molecule or a nanowire branch leads both to a reduction, Ak, of the 
thermal conductance and a reduction, AG, of the electronic conductance. As we show in Ref. 
|12j . at room temperature (RT) the ratio An/ AG > 50 for the alkyl functionalized SiNWs, and 
An/ AG > 20 for a nanotree. By engineering the SiNW surfaces it is thus possible to reduce the 
phonon conductance while keeping the electronic conductance almost unaffected. Such surface 
decorated SiNWs would be promising candidates for nanoscale thermoelectric applications. 
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Figure 1. (Top) Sketch of electronic density of states in the alkyl functionalized SiNW (middle) 
and in the nanotree (bottom). The trunk in the nanotree is broader than the branch and thus 
has a smaller bandgap. For both structures electrons and holes close to the band edges are 
weakly scattered while phonons are strongly scattered. 



2.3. Systems 

We consider two specific systems shown in Fig. [TJ The first is an alkyl functionalized SiNW 
with a wire diameter of 12 A, and with the wire oriented along the (110) direction. The alkyl 
(C n H2 n +i) is attached to the H-passivated nanowire replacing a H atom. The second system 
is a nanotree, where a small diameter (12 A) branch is attached to a larger diameter (20 A) 
trunk. The trunk is oriented in the (110) direction while the branch is oriented along the (110) 
direction, and is thus perpendicular to the trunk. The length of the branch, Lb, is varied. 



2.4- Methods 

The electronic Hamiltonian, H, and overlap matrix, S, of the alkyl functionalized SiNWs are 
obtained from local orbital DFT calculations [22]. The calculations are performed on super-cells 
containing 7 wire unit cells, as indicated in Fig. [TJ 

For the nanotrees, we use a tight-binding (TB) model since these systems contain > 
1100 atoms, too many for our DFT implementation. The electronic TB Hamiltonian 
describing the nanotree is calculated using a 10 band sp s d 5 s* nearest-neighbor orthogonal TB 
parametrization [23, 24J. We recently applied the same TB model to study thermoelectric 
properties of surface disordered SiNWs [7]. The same TB parameters were recently also applied 
to study SiNW band structures [25J and surface roughness [10|, 126], 

The phononic system is described using the Tersoff empirical potential (TEP) model [271 128] 
for both the nanotree and the functionalized SiNW. For pristine wires, we have recently shown 
that the TEP model agrees well with more elaborate DFT calculations [11]. We limit our 
description to the harmonic approximation, thus neglecting phonon-phonon scattering. The 
harmonic approximation is always valid at low temperatures. In bulk Si, the room temperature 
anharmonic phonon-phonon relaxation length at the highest frequencies is \ a {u max ) ~ 20 nm 
and increases as A a oc u~ 2 at lower frequencies |29] , Experimental studies of silicon films |30] 
showed that the effective mean free path of the dominant phonons at room temperature is ~300 
nm. For relatively short wires with lengths L < 100 nm the anharmonic effects thus seem to be 
of limited importance, and the harmonic approximation is expected to be good. 

We calculate the electronic conductance from the electronic transmission function, T e (s) 
following the standard non-equilibrium Green's function (NEGF)/Landauer setup, where the 
scattering region (i.e. the regions shown in Fig. [1]) is coupled to semi-infinite, perfect wires |31j . 
Explicitly, the transmission function is given by 

T e (e) = Tr[G r (e)T L (e)G a (e)T R (e)}, (1) 

where the retarded and advanced Green functions G r ' a can be calculated recursively, and the 
coupling matrices Tl^r are evaluated with a DFT calculation. A discussion of the validity of this 
expression can be found, e.g., in Ref.|31j. The electronic quantities in the ZT formula can then 
be written as [321E3JIM] G e = e 2 L , S = Li(/i)/[eTL (^)] and n e = [L 2 (p)- (^i(^)) 2 / L (jj)]/T 
where L m (n) is given by 

LM = I £L d£%{£){£ - » )m (-^fe^) • (2) 

The phonon transmission function, 7^/i(w), at frequency cv is calculated in a similar way as 
the electronic transmission with the substitutions H — > K (K is the force constant matrix) and 
eS —J- w 2 M, where M is a diagonal matrix with the atomic masses [35} l36| 137] . The phonon 
thermal conductance is 

7j2 poo e hoj/k B T 

^ (T) = 2^t4 d ^ 2 W (e ^ B T_ 1)2 " (3) 

2.5. Charge transport 

Figures [2] (a) and (b) show the calculated hole and electron transmissions for a pentyl (C5H11) 
functionalized SiNW. Notably, the transmission is nearly perfect close to the bandedges, in 
agreement with the findings of Ref. [2]. Figures 2 (c) and (d) show the projected density of 
states (PDOS) on the wire and on the pentyl. The high transmission regions in panel (a) and 
(b) are seen to correspond with regions of vanishing PDOS on the pentyl molecule. Likewise, 





Figure 2. Hole and electron transmissions for pentyl functionalized, D = 12 A SiNW (a)-(b) 
and for a nanotree with D = 20 A trunk and D = 12 A branch (e)-(f). Panels (c)-(d) and (g)-(h) 
show the PDOS on the wire and on the pentyl and branch, respectively. The energy scales are 
relative to the valence band edge, E v for the holes (left column) and relative to the conduction 
band edge, E c for electrons (right column). The DFT band gap of the D = 12 A wire is 1.65 eV, 
while the TB band gap of the D = 20 A wire is 1.77 eV (from Ref.[]2]). 



at energies in the valence band where scattering is observed, there is a relatively large PDOS at 
the pentyl. 

Figures 2 (e) and (f) show the transmission through the nanotree. The branch length is 
15.4 A. Again, the transmission close to the band edges is nearly perfect. Figures 2 (g) and 



(h) show the PDOS on the main wire and on the branch. We again observe a correspondence 
between perfect transmission and low PDOS on the branch. 

The almost perfect transmissions close to the band edges can be qualitatively understood 
from the schematic drawing in Fig. [1] (top). The HOMO and LUMO level of the pentyl are 
located deep inside the bands |14| and the molecular states are thus not accessible for electrons 
or holes close to the band edges. For the nanotree, the branch has a smaller diameter and thus 
a larger bandgap. Electrons or holes in the trunk, with energies close to the band edges, are 
not energetically allowed in the branch and therefore do not 'see' the branch. In addition to 
the energy considerations, the spatial distribution of the Bloch state also plays a role: the first 
valence and conduction band Bloch states of the main wire have more weight in the center of 
the wire than at the edge [7] . 

2.6. Thermal transport 

Figure [3] shows the temperature dependence of the thermal conductance ratios k/kq, where kq is 
the pristine wire thermal conductance, which in the low energy limit equals the universal thermal 
conductance quantum, kq(T) = 4(ir 2 k B T/3h) [38]. Panel (a) shows the ratios for wires with 
alkyles, C n H2„+i, with different lengths, n = 3,5, 7. The thermal conductance at RT is reduced 
by ~ 10%, and the overall behavior does not depend on the alkyl length. The inset shows 
the phonon transmission at low phonon energies. Note the resonant dips in the transmission, 
where exactly one channel is closed yielding a transmission of three. These dips are associated 
with an increased local phonon density of states at the alkyl molecule at the resonant energies, 
corresponding to a localized vibrational mode. Such Fano-like resonant scattering is well-known 
from electron transport [39] . A phonon eigenchannel analysis [3D] shows that the transmission 
dips are due to a complete blocking of the rotational mode in the wire. The corresponding 
localized alkyl phonon mode is a vibration in the plane perpendicular to the wire axis. 

Panel (b) in Fig. shows the thermal conductance ratio for nanotrees with different branch 
length, Lb- There is only a weak dependence on Lb at low temperatures, and at RT the 
four curves basically coincide showing a thermal conductance reduction of 17% of the nanotree 
compared to the pristine wire. Again we observe resonant transmission dips for the nanotrees. 
Two channels - the rotational and one flexural mode - close completely at the resonance due to 
two quasi-localized vibrational modes in the branch. These phonon backscattering resonances 
are responsible for the dip in the k/kq ratio around T = 10 K. Notice that all the conductance 
ratios approach unity in the low temperature limit. This is because the four acoustic modes 
transmit perfectly in the limit u — > [38j. 

We may vary the thermoelectric figure of merit, ZT, by varying the chemical potential. 
Typically ZT displays a maximum for \i close to the band edge [7J [5] . Figure 0] shows the 
maximum ZT values for the pentyl functionalized SiNW (squares), the nanotree (circles), and 
surface disordered SiNWs (triangles), where disorder is modeled by introducing surface silicon 
vacancies. The diameter of the surface disordered wire is D = 20 A and it is oriented along 
the (110) direction. The calculational details are given in Ref. [JJ. The curves show ZT as a 
function of the number (N) of pentyl molecules/nanotree branches/silicon vacancies. 

Figure [3] shows that increasing the number of scattering centers, i.e. the number of pentyl 
molecules or nanotree branches increases the ZT for both hole transport (a) and electron 
transport (b). In the case of holes in the pentyl functionalized SiNWs, the ZT reaches an 
almost constant level of ZT = 0.7 at N = 40, but in all other cases ZT increases throughout 
the range. Increasing the density of molecules/nanotree branches or increasing the length of the 
wire will thus increase the thermoelectric performance. The surface disordered wires (triangles) 
show an increasing ZT vs N but at values significantly lower than the two other surface modified 
wires. 
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Figure 3. Thermal conductance ratio k/kq for (a) pentyl functionalized SiNWs with different 
alkyl lengths, and (b) for nanotrees with different branch lengths, Lb- The nanotree trunks are 
again D = 20 A with 12 A diameter branches. The insets show the phonon transmission function 
at low energies. Fano-like resonant scattering is observed in both systems (from Rcf. [12 ). 



2.7. Discussion 

A number of idealizations have been made in our calculations, and we next assess their 
significance. The structures we have considered represent plausible choices, dictated by 
computational limitations, but do not necessarily match quantitatively real structures. Thus, 
for example, a surface decorated SiNW will also be rough, and one should consider the combined 
effect of all scattering mechanisms. We have not carried out optimizations neither with respect 
to the attached molecules nor with respect to the geometry of the nanotrees. Electron-phonon 
and phonon-phonon scattering will affect both the electronic and thermal conductances and the 
obtained ZT values [SJ. We do not expect to reach quantitative agreement with experiment but 
believe to have identified important trends: In SiNW based thermoelectrics, surface decorations 
in terms of added molecules or nanowire branches seem to be a better approach than surface 
disorder in the ultra-thin limit. 




Figure 4. Thermoelectric figure of merit, ZT, for p-type (a) and n-type ( b) wires. N is the 
number of pentyl molecules (squares), nanotree branches (circles) and silicon surface vacancies 
(triangles) in the wire (from Ref. [12J). 



3. Microscopic theory for damping of vibrations in atomic gold wires 

3.1. Introductory remarks 

As pointed out in the beginning of this review, a particularly important issue in molecular 
electronics concerns the energy exchange between the charge carriers and the molecular contact. 
Thus, the local Joule heating resulting from the current passing through the contact, and its 
implications to the structural stability of such contacts are presently under intense investigation 
|41|. 142 1 l%3"l 1441 145| . Experimentally, local heating in molecular conductors in the presence of the 
current has been inferred using two- level fluctuations [46J and Raman spectroscopy [47 . 

Mono-atomic chains of metal atoms |48| are among the simplest possible atomic-scale 
conductors. The atomic gold chain is probably the best studied atomic-sized conductor, and a 
great deal of detailed information is available from experiments and related theoretical studies; 
see Ref. [51] for an extensive list of references. 

In the case of a gold chain Agrait et al [50] reported well-defined inelastic signals in the 
current-voltage characteristics. These signals were seen as a sharp 1% drop of the conductance 
at the on-set of back-scattering due to vibrational excitation when the voltage equals the 
vibrational energy. Especially for the longer chains (6-7 atoms), the vibrational signal due 
to the Alternating Bond-Length (ABL) mode [52, 55J, dominates. This resembles the situation 



of an infinite chain with a half-filled electronic band where only the zone-boundary phonon can 
back-scatter electrons [19] due to momentum conservation. 

The inelastic signal gives a direct insight into how the frequency of the ABL-mode depends 
on the strain of the atomic chain. This frequency can also be used to infer the bond strength. 
The signature of heating of the vibrational mode is the non-zero slope of the conductance versus 
voltage beyond the on-set of excitation: with no heating the curve would be flat. Fits to the 
experiment on gold chains using a simple model [S3] suggest that the damping of the excitation 
can be significant. However, the experiments in general show a variety of behaviors and it is not 
easy to infer the extent of localization of the ABL vibration or its damping in these systems^. 

In order to address the steady-state effective temperature of the biased atomic gold chain 
theoretically, it is necessary to consider the various damping mechanisms affecting the localized 
vibrations, such as their coupling to the vibrations in the contact, or to the phonons in the 
surrounding bulk reservoirs. Here we outline how to calculate the vibrational modes in atomic 
gold chains and their coupling and the resulting damping due to the phonon system in the leads; 
a full account is available in our recent paper [59] . We work within the harmonic approximation 
and employ first principles density functional theory (DFT) for the atomic chain and the contacts 
[60] while a potential model is used for the force constants of the leads |61| . 

Here we focus on chains between two (lOO)-surfaces. We consider chain-lengths of 3-7 atoms 
and study the behavior of their vibrations and damping when the chains are stretched. TEM 
micrographs indicate that the chains are suspended between pyramids, so in our calculations we 
add the smallest possible fcc-stacked pyramid to link the chain to the given surfaces. 

3.2. Method 

In the microscopic modeling of the dynamics of the atoms forming the junction one must 
distinguish between different kinds of atoms: those belonging to the chain, those belonging 
to the pyramid joining the chain to the substrate, and substrate atoms. The technical issues of 
how these different types of atoms are treated are discussed in Ref.[59]. A quantity of central 
importance to all analysis is the mass-scaled dynamical matrix, K, 



where E is the total energy of the system, Uj is the coordinate corresponding to the z'th degree 
of translational freedom for the atoms of the system. rrii is the mass of the atom that the z'th 
degree freedom belongs to. K governs the the evolution of the vibrational system within the 
harmonic approximation. The Newton equation of motion reads 



where A denotes a mode of oscillation in the system and e\ is the corresponding quantization 
energy. The evaluation of K proceeds as follows. Finite difference DFT calculations were used 
for the nanostructure, while for the surfaces we used an empirical model due to Treglia and 
Desjonqures [61J. The position of the interface between the region treated by DFT and the 
region treated by the empirical model is a parameter that can be varied; here one must find a 
value for which the physical results converge keeping the computational labor manageable. 

3. 3. Green's functions projected on the atomic chain 

All properties of interest in the present context can be derived from the Green's function D, 
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1 N. Agrai't, private communication. 



where rj = + and we denned the inverse of the Green's function by M = D _1 . Specifically, we 
shall need the Green's function T)dd projected onto the region close to the atomic chain. Our 
procedure is based on a method due to Mingo et al [62] which has previously been tested in an 
investigation of finite Si nanowires between Si surfaces. We find [59] 

Unperturbed 




Figure 5. [Top] The forces between surface atoms within next-nearest neighbor distance (4. 08 A) 
of the added atoms are perturbed by the presence of the added atoms. [Bottom] The device 
region is where the coupling between the atoms is different from the values for the two 
unperturbed surfaces. The coupling between the device region and the leads is considered 
to be unperturbed. 



T) DD = [M DD -M Da (M aQ )" 1 M QD ]- 1 

= [Mdd-Udd]- 1 , (7) 

where the index a = (L,R), i.e., the left and right unperturbed surface, while D = {1,A, 2} 
labels the regions in left contact, the atomic chain, and the right contact, that are different from 
the bulk (see Fig. This expression defines the self-energy IId£> = M.£) a (M. aa )~ 1 'M. a £) which 
has the matrix structure 

/ M 1L (M LL )- 1 M L1 \ 

U DD = . (8) 

V M 2R (M R fi)- 1 Mfi 2 / 

This object can be evaluated as follows. First, in the limit of large regions 1 and 2, the coupling 
elements M^i and M# 2 must approach those of the unperturbed surface, Mf^ and M^ 2) 
respectively. In what follows, we shall make the approximation that the regions 1 and 2 are 
chosen so, that this condition is satisfied sufficiently accurately. Second, we note that the 
matrix M aa is indistinguishable from the matrix Mf a , as long as the involved atoms are outside 
the perturbed regions 1 or 2. Therefore, we can write 

M 1L (M LL )- 1 M L1 ~ Mf L (Mf L )- 1 Mf 1 = nfi 
M 2R (M RR )- 1 M R2 ~ Mf R (M| JJ )- 1 Mg 2 = nf 2 , 

(9) 

where the accuracy increases with increasing size of regions 1 and 2. On the other hand, using 
the definition of the self-energy, we have 

nfx = Mf 1 -(Df 1 )- 1 

nf 2 = Mf 2 -(Df 2 )-\ (10) 



Figure 6. Distances and angles used to define the average bond length, B = (bj), and the 
average bond angle, T = (6j), respectively 
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Figure 7. Force as function of average bond length B. 



where D^, i = 1, 2 is the projection of the unperturbed Green's functions onto the atoms in regions 
1,2, respectively. This object is evaluated by exploiting the periodicity in the ideal surface plane. 
The Fourier transform of M s in the parallel directions has a tridiagonal block structure and we 
can solve for its inverse very effectively using recursive techniques (see e.g. Sancho et al [63J). 
An analysis of the convergence properties of this procedure is given in Ref. [59] . 

To sum up, the calculation is preformed in the following steps: (i) Start with perfect leads 
and specify the device in between them, (ii) The atoms in the leads where K is perturbed by the 
presence of the device are identified, (iii) The unperturbed surface Green's function D 5 is found 
via fc-point sampling and then used to construct the self-energy, Eqs.([9T H0]) . (iv) The perturbed 
Green's function is then found using this self-energy via Eqs.flJHHJ). With the self-energy at 
hand, we can define the life-times of the various modes in a standard fashion (see Ref. [59J for 
details), and we now proceed to present some of our results. 

4. Results 

4-1- Geometrical Structure, the Dynamical Matrix, and Q-f actors 

Here we report results for the geometrical structure of the chains as a function of the chain 
elongation. For each geometry (identified by the number of atoms in the chain, the surface 
orientation and the type of contact to substrate) a range of calculations were performed with 
the two surfaces at different separations, with the separations incremented in equally spaced 
steps. To be able to compare chains of different lengths and between different surfaces we define 
the average bond length, B = (bj), as the average length between neighboring atoms within the 
chain, where j runs over the number of bonds in the chain (see Fig. [6]). B is useful because it is 
closely related to the experimentally measurable force [51] on the chain. The close relationship 
between B and the force is demonstrated in Fig. We note that the force vs. B curves to a good 
approximation follows a straight line with a slope of k = 2.5 eV/A , which can be interpreted 
as the spring constant of the bonds in the chain. 

As the systems are stretched it is mostly the bonds in the chain that are enlongated. At low 
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Figure 8. Average bond angle T = (8j) as a function of the average bond length. The long 
chains adopt a zig-zag structure (T < 180°) at low B while the short chains remain linear. 
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Figure 9. Local energies of a 4 atom chain between two (lOO)-surfaces at different strains. The 
largest eigenvalues are connected by a line to guide the eye. 



B Fig. [8] shows that the longer chains adopt a zig-zag confirmation at low average bond length. 
The 3- and 4-atom chains, however, remain linear within the investigated range. Furthermore, 
the longer chains have a similar variation in the average bond angle. Our results justify the 
earlier theoretical studies by Frederiksen et al [33] and Sanchez-Portal et al |64| . 

We next investigate the energies that are related to different types of movement by analyzing 
the eigenmodes and eigenvalues of selected blocks of K. In particular, we can consider the local 
motion of individual atoms or groups of atoms, freezing all other degrees of freedom, by picking 
the corresponding parts of K. The square root of the positive eigenvalues of the reduced matrix, 
which we call local energies, give the approximate energy of a solution to the full K that has 
a large overlap with the corresponding eigenmode, if the coupling to the rest of the dynamical 
matrix is low. The negative eigenvalues of a block are ignored since they correspond to motion 
that is only stabilized by degrees of freedom outside the block. 

The behavior of the dynamical matrix in terms of local energies is illustrated in Fig. [9) In the 
chain the local energies are quickly reduced with increased strain while the dynamical matrix 
of the surfaces is hardly affected at all. The middle bonds in the chain are the ones that are 
strained and weakened most when the surfaces are moved apart. It is also where the chain is 
expected to break [53] . 

Fig. [TO] shows how motion parallel to the chain occurs at higher energies than perpendicular 
motion, and that the LO type motion of the ABL modes has the highest energy. The local 
energies of the ABL/LO type motion move past the local energies of the pyramid as the strain 
is increased. In this way the ABL/LO modes can in some sense act as a probe of the contacts. 

Figure Fig. [TT] shows the projected DOS for a chain with 4 atoms at an intermediate strain. 
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Figure 10. Local energies for selected blocks of K as a function of the average bond length 
in the 4 atom chain. Since the central chain in this case consists of 2 atoms we can classify 
the eigenvectors as LO: longitudinal optical, LA: longitudinal acoustic, TO: transverse optical 
(doubly degenerate) or TA: transverse acoustic (doubly degenerate). 



Notice the large variation in the width of the peaks. The peaks with a small width correspond 
to modes that have the largest amplitude in the chain, while the peaks with a large width 
correspond to modes with large amplitude on the interface. 




Figure 11. Projected DOS onto a representative selection of the vibrational modes of the device 
region. 

In Fig.[T2]we present the Q-factor, spatial localization and peak energy of all modes for chains 
with 3-7 atoms between (100) surfaces. The ABL modes are of special interest. These modes 
have been identified by previous theoretical and experimental studies as the primary scatterers 
of electrons jBl H EZ1 I53" l 157 ] . The ABL modes are easily identified in Fig. [T2l since they 
have the highest energy of the modes that are spatially localized to the central chain (black or 
dark gray on the figure). Modes corresponding to transverse motion of the central chain are also 
clearly visible. These modes are energetically and spatially localized, but are of limited interest 
because of a low electron-phonon coupling. 

Certain ABL modes are very long-lived. At low strains, ABL modes lie outside the bulk 
(and surface) band and have, in the harmonic approximation, an infinite Q-factor. In reality the 
Q-factor will be limited by electron-phonon and anharmonic interactions. At higher strain the 
ABL modes move inside the bulk band and one observes a great variation in the corresponding 
Q-factors. When the peak energy lies inside the bulk band it will mostly be the structure of the 
connection between the bulk crystal and the chain that determines the width of the peak. 
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Figure 12. The vibrational modes for chains with 3-7 atoms between two 100-surfaces. The 
center of the disks are positioned at the peak of the projection of vibrational DOS on the mode 
in question. The area of a disk is proportional to the Q\, but is limited to what corresponds to 
a Q-factor of 250. The gray level is a measure of the localization of the mode, black being most 
localized. 



The long chains tend to have longer lived ABL /LO modes due to the larger ratio between the 
size of the chain and the size of its boundary. The 7-atom chain is especially interesting since it 
has an ABL/LO type mode with a damping of 5 meV at one strain, while at another strain the 
ABL/LO mode has a damping of 300 meV, i.e., more than one order-of-magnitude variation in 
the harmonic damping of the primary scatterer of electrons due to only a 0.03 A change in the 
average bond length! 

The largest damping of an ABL-mode for these systems is 7a ~ 1 meV, which is still 
significantly lower than the ~ 20 meV band width. This can be attributed to fact, noted 
above, that there always exists a large mismatch in local energies moving from the central part 
of the chain to the rest of the system (see Fig. [U]) . 

Previous studies by Frederiksen et al [55J obtained a rough estimate for the variation of the 
non-electronic (harmonic and anharmonic) damping of 5-50 /ueV for the longer chains by fitting 
the experimental IETS signals of Agra'it et al [IS] to a model calculation. The estimated peak 
energies lie well within the bulk band for all the recorded signals. The reason the non-electronic 
damping rate can be extracted is because the excitation of vibrations and damping of vibrations 
through electron-hole creation are both proportional to the strength of the electron-phonon 



coupling. This means that the step in the experimental conductance, when the bias reaches the 
phonon energy, can be used to estimate strength of the electron-phonon interaction and thereby 
the electron-hole pair damping. The slope in the conductance beyond this step can then be 
used to extract the total damping. By subtracting the electron-hole pair damping from the 
total damping we get an estimate of the sum of harmonic and anharmonic contributions to the 
damping. 

The estimate in Ref. [SS] agrees well with our lowest damping of 5 //eV. The highest damping 
we have found was ~ 400 /xeV found for the 6 atom chain which is an order of magnitude larger 
than the upper limit of Ref. [55J. We believe that this discrepancy can be largely attributed 
to the difficulty in extracting the necessary parameters from experiments when the harmonic 
damping is large. Furthermore, for the 6-7 atom chains we observe that the high damping occurs 
at low strain, where the electron-phonon coupling is weak |49| . 



5. Conclusion and Discussion 

Our most important finding is that the values of the harmonic damping for the ABL modes can 
vary by over an order of magnitude with strain. Even with small variations in the strain, the 
harmonic damping can exhibit this strong variation. This extreme sensitivity may explain the 
large variations seen experimentally in different chains. 

The range of the harmonic damping also depends strongly on the number of atoms in the 
chain since we see a clear increase in localization going from a 6- to a 7-atom chain. The chain 
with 7 atoms really stands out, since it, in addition to having very localized modes in general, 
it also has the greatest variation in harmonic damping. This strong variation in the harmonic 
damping of the ABL /LO-modes suggests that accurate atomistic calculations of the vibrational 
structure is necessary to predict the inelastic signal. 

The techniques described here can be combined with electronic transport calculations to 
predict the inelastic signal in the IV characteristic of a system. This will be done in future 
work, where we will also eliminate the use of the empirical model for the leads and use ab-initio 
parameters for the entire system. 
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